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(Nl ■ ABSTRACT 

We show that in a system of two planets initiahy in nearly circular orbits, an impulse 
perturbation that imparts a finite eccentricity to one planet's orbit causes the other 
planet's orbit to become eccentric as well, and also naturally results in a libration of 
, their relative apsidal longitudes for a wide range of initial conditions. We suggest that 

^Sj I such a mechanism may explain orbital eccentricities and apsidal resonance in some 

^ exo-planetary systems. The eccentricity impulse could be caused by the ejection of a 

Q ' planet from these systems, or by torques from a primordial gas disk. The amplitude of 

■ secular variations provides an observational constraint on the dynamical history of such 

systems. 

Q . Subject headings: planetary systems — celestial mechanics — instabilities - stars:individual(f 

^ ; Andromedae, HD 83443) 



Introduction 



X 

. 5t , Orbital eccentricities in exo-planetary systems^ discovered thus far are often surprisingly large 

and have proven to be a major puzzle in understanding these systems. At least two systems with 
multiple planets in eccentric orbits {v Andromedae and HD 83443) are suspected of exhibiting 
secular apsidal resonance (Chiang, Tabachnik & Tremaine 2001, Wu & Goldreich 2002). Apsidal 
resonance is the phenomenon of phase-locking of the apsidal longitudes of two orbits, such that the 
two planets have a common average rate of apsidal precession and the angular difference of their 
apsidal longitudes, Aw = wi — W2, librates around 0. 

We describe here a dynamical mechanism for establishing apsidal resonance in a pair of planets 
that are initially on nearly circular orbits. We use classical analysis from celestial mechanics to show 
that an impulse perturbation that imparts an eccentricity to one of the orbits, excites the other 
planet's eccentricity on a secular timescale and also results in the libration of the relative apsidal 



^http: / /www. exoplanets.org 
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longitude for a wide range of initial conditions. A plausible cause of an impulse perturbation is the 
ejection of a planet from the system; torques from an exterior primordial disk may also cause such 
an eccentricity impulse. 



2. Secular dynamics 

Consider a pair of planets in well-separated orbits, whose orbital period ratio is not too close 
to a ratio of two small integers. Under the assumption that the relative inclination of their orbital 
planes is small, the secular interactions in such a system preserve the orbital semimajor axes and 
induce quasi-periodic variations in the other orbital parameters. For small eccentricities and small 
relative inclinations, the ccccntricity-apsidal line variations arc decoupled from the inclination- 
nodal line variations (Murray & Dcrmott 1999). Under these conditions, the secular variations 
of the eccentricities and apsidal longitudes are described by a set of first order linear differential 
equations for the components of the so-called eccentricity vectors, {kj {t) , hj (t) ) = ej (cos Wj , sin xuj ) , 
for each of the planets, where ej is the eccentricity and Wj is the apsidal longitude. The general 
solution for two planets is a superposition of two eigenmodes, 

hjit) = Efhm{git + I3i)+Efs{n{g2t + P2), 

kj{t) = Efcos{git + pi) + Efcos{g2t + /32). (1) 

Let myt, mi, m2 be the masses of the star, the inner planet and the outer planet, respectively; ni, 
n2 the orbital mean motions of the inner and outer planet, respectively; a = ai/a2 < 1 the ratio of 
the semimajor axes, bs\a) a Laplace coefficient. We also define jj, = mi/m2, and B = ^3/2/^3/2" 
With these definitions, the mode frequencies are 

51,2 = ^\/a^3/2(a) J[J^«2{l + /xv^ ± [(1 - fiV^f + 4/iv^s2] ^^^|, (2) 

and the eigenvector amplitudes E^''^ are given, up to a normalization factor, by 

E^^^^ _ 1 - T [(1 - /^V^)' + 4/xVSS2]i/2 



^(i) 2B 



where the '— ' sign in ^ corresponds to j = 1, and the sign to j = 2. The phases Pj and the 
normalization factors for the amplitudes of the modes are determined by initial conditions. 

We assume initially nearly circular orbits, as is expected of planet formation in a dissipative 
protoplanetary disk. And we consider an impulse perturbation^ that imparts a finite eccentricity. 



^In general, an impulse perturbation would result in not only an eccentricity impulse but also impulse perturbation 
to other orbital elements. The secular dynamics described here remains valid provided we use the post-impulse value 
of a2 in evaluating a and n2 ■ 
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e2f , to the outer planet's orbit, at time t = 0. (The case of an impulse to the inner planet's orbit is 
solved in a similar way.) Then, with initial conditions hi{0) = /iio,fei(0) = A;io,/i2(0) = 0, A;2(0) = 
e2f, we have 

77.(1) P'^^W - e2f „(1) . a P2hw 

E{ ' cos /?i = , E\ ' sm/3i — 



Pi - Pi P2- Pi 

17.(2) a Plkl0-e2f „(2) . a Pi ^10 /.N 

E\'cos(32 = , E\'sml32 = • (4) 

Pi- P2 Pi- P2 

This solution admits either circulation or libration of the relative apsidal longitude, Aw. Examples 

of three types of possible behavior are shown in Fig. 1. In each of these cases, the cumulative 
probability density function of |Atu| and of ei (i.e., the fraction of time that |Atu|,ei is less than 
some value) is also shown. 

In the case of zero initial eccentricity, eio = [Fig. l(a-d)], mode 1 and mode 2 appear with 
equal amplitude but opposite phase in the inner planet's eccentricity vector. The time evolution of 
the planets' eccentricities and of Aw is given by 

1/2 



ei{t) = I ^^^^ I 1 - cos(5i - g2)t 
\Pi - P2 



e2{t) = : — — — 7 p\+ pI- 2piP2 cos(5i - g2)t 

1/31 -P2|L J 

imiAwit) = ^^^tan[-(5i-52)t + J]. (5) 

P2+ Pi ^ 

The eccentricity of the inner planet grows from its initial value of zero to a value comparable to 
e2f on the secular timescale ~ l^i — 52!"^; the eccentricities vary periodically (but not sinusoidally) 
with a common frequency, gi — g2, and with opposite phases. The inner planet's eccentricity has 
a maximum excursion to 2e2f/|/5i — P2\, while the eccentricity of the outer planet varies between a 
minimum of e2f |pi + /92|/|pi — P2I and a maximum of e2f . The apsidal longitude difference is limited 
to (—90°, 90°); it oscillates with the same frequency, gi — g2- The singularity of the tangent function 
at i^TT implies a discontinuity in the time evolution of Aw. This discontinuity is not unphysical 
(see Fig. 2 below): as the eccentricity vector of the inner planet passes periodically through the 
origin, its apsidal longitude changes discontinuously. 

A small non-zero initial eccentricity, eio <C e2f, results in apsidal libration with an amplitude 
smaller than 90° for values of wiq that are approximately within ±90° of W20 [e.g., Fig. l(e-h)], 
and apsidal circulation otherwise [e.g.. Fig. l(i-l)]. In all cases, the secular dynamics following an 
impulse perturbation to 62 yields a low probability of finding \Aw\ > 90° [see Fig. l(c,g,k)], and 
only ~ 1/3 probability of finding e\ less than half its maximum [see Fig. l(d,h,l)]. As we show 
below, a small non-zero eio results in apsidal resonance with approximately 50% probability. 

It is illustrative to consider the limiting case m\ — > 0. In this case, the orbital elements of 
the outer planet are constant: k2{t) = e2f,/i2(i) = 0; for the inner planet, we have the following 
solution: 



1/2 



ki{t) = +{kio — eif)cosgt + hiosmgt + eif, 
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hi (t) = - {kio - en) sin gt + hw cos gt, (6) 

with ^ 

eif = Be2{ and g = -ni — -a^b^^Lla). (7) 

4 

This solution admits a simple geometrical interpretation as illustrated in Fig. 2a: the inner planet's 
eccentricity is the vector sum of a 'forced' component of magnitude eu which is aligned with the 
outer planet's apsidal longitude (taken to be zero here, without loss of generality), and a 'free' 
component of magnitude 

ei,free = (e?o + e?f - 2eif /cio)^''^ (8) 

which rotates with frequency g. Note that eif is independent of the mass of the outer planet, but 
the frequency g is proportional to that mass. 

Apsidal libration occurs if ei^free < eif. In the (A;io,^io) plane, this condition is satisfied in 
the area of intersection of two circles, one of radius eio centered at the origin, the other of radius 

eif centered at (eif,0), as shown in Fig. 2b. Assuming random values of zuiq and eio < £, the 
probability P^b that the impulse on the outer planet will result in apsidal libration is simply the 
intersection area as a fraction of the area of the circle of radius e. This is readily evaluated: 

For e <C £21, we have Piib = 0.5: nearly all values of wiq between —90° and +90° result in apsidal 
resonance. For e = en, the probability of apsidal resonance drops only slightly, to ~0.4. 

For the general case, mi 7^ 0, the probability of apsidal resonance can be determined numer- 
ically. It is easy to see from the general solution (eq. 4) that Pub has only a weak dependence on 
the mass ratio fi = mi/m2, and Pub remains close to 0.5. 

When apsidal libration occurs, the libration amplitude is a steep function of zuiq and eio/e2f. 
This is illustrated in Fig. 3, for one choice of parameters, /j, = 0.5 and a = 0.33, corresponding to 
the case of planets C and D of v Andromedae. 



3. Discussion 

We have shown that in a two planet system, an impulse perturbation on the eccentricity of the 
outer planet can excite the eccentricity of the inner planet on a secular timescale, and will result in 
apsidal libration with a ~ 50% probability. There are two significant characteristics of the secular 
dynamics following an eccentricity impulse, whether or not it results in apsidal resonance: (i) the 
inner planet's eccentricity has large amplitude variation, and it drops to its initial small value 
periodically; (ii) only a small fraction of time is spent at large values of |Atu| and small values of 
ei (see Fig. 1). Therefore, observing small values of |Atu| and large values of ei is not surprising if 
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the system has suffered an impulse perturbation in its history and the secular dynamics indicates 
large amplitude variations. 

We note that in the opposite limit of a slow adiabatic perturbation which increases one planet's 
eccentricity on a timcscalc much longer than the secular timescale, apsidal resonance will occur with 
nearly 100% probability and the resulting libration amplitude will be very small (for initially circular 
orbits) . It is rather remarkable that the secular dynamics of two planets admits a high probability 
of apsidal resonance in both the impulse and the adiabatic limit of eccentricity perturbation of one 
planet. We defer detailed analysis to a future study, but we note here that the apsidal libration 
amplitude provides an observational diagnostic: large amplitudes favor the impulse mechanism 
whereas small amplitudes favor the adiabatic mechanism. 

One possible mechanism for the eccentricity impulse is planet-planet scattering. Highly eccen- 
tric orbits can be produced by gravitational interactions of planets initially on circular orbits that 
are close to the threshold of orbital instability. This has been proposed as a possible explanation 
for the eccentric orbits of exo-planets (Rasio & Ford 1996, Weidenschilling & Marzari 1996, Lin &: 
Ida 1997, Ford et al. 2001, Marzari &; Weidenschilling 2002). There are three possible outcomes of 
such interactions: a final stable configuration different from initial, the merger of two planets, or 
the ejection of one planet. Numerical experiments thus far have explored only a small fraction of 
the parameter space of two and three giant-planet systems. In systems of two equal mass planets, 
comparable branching ratios for the three possible outcomes are found; the two former outcomes 
do not result in eccentricity excitation, but the third - ejection of a planet - does, with typical 
eccentricity of the surviving planet in the range 0.4-0.8 (Ford et al. 2001). One problematic aspect 
is that two planet systems must have initial orbits fine-tuned close to the threshold of dynamical 
instability, i.e. orbital separation close to 2\/3 times their mutual Hill radius (Gladman 1993). But 
this problem does not exist in systems of three or more planets, where chaotic instabilities can 
occur over a wider range of orbital separations; such systems can persist in quasi-stable orbital 
configurations for long periods of time, possibly exceeding 10^ years, before becoming dynamically 
unstable; then the most common outcome is the ejection of one planet (Marzari & Weidenschilling 
2002). Also, in the more realistic case of unequal masses, the probability of planet ejection is likely 
to be enhanced further. 

We consider the ejection into an unbound orbit of a third planet by gravitational scattering from 
the outer of the two planets whose secular dynamics we have analyzed in section 2. (For comparable 

masses, ejection by scattering from the outer of the two planets is more likely than from the inner 
planet because the escape velocity is smaller at larger distance from the star: Vqsc ~ r~^/^; it is also 
likely to produce the least direct perturbation to the inner members of the system.) Assuming that 
the ejected planet leaves on a nearly parabolic orbit, energy and angular momentum conservation 
yield the following estimate for the final eccentricity of the bound planet: 



where 7713,03,53 are the mass, initial semimajor axis and final periastron distance of the ejected 




(10) 
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planet; is the initial semimajor axis of m2, and is related to its final semimajor axis by energy- 
conservation: ~ 'm-20'2(i3 / + "1-3^2) • Unfortunately, there is no additional dynamical 
constraint or simple argument to constrain the final periastron distance of the ejected planet. So 
the final eccentricity of the bound planet cannot be predicted analytically and must be determined 
numerically. Still, we can make the following estimate: since the periastron distance of the ejected 
planet is likely to be not greatly in excess of the apoastron distance of the surviving planet, then, 
for §3/02 ^ 1-1.5, the final eccentricity of the bound planet will be in the range of 0.2-0.7 for 
m^/m2 — 0.1-0.7. 

Once a dynamical instability sets in, the subsequent evolution is highly chaotic and unpre- 
dictable in detail. A simplified description, based upon the presently available numerical simula- 
tions, is as follows: at first, close encounters between m2 and 771,3 Isad to a rapid increase in their 
orbital eccentricities, over ~ 10 — 10^ orbital periods; this is followed by a longer period, > 10^ 
orbits, during which many weak distant encounters gradually and stochastically increase the orbital 
period and apoastron of m^, the planet-to-be-ejected. In the present context, the initial eccentricity 
excitation of the surviving planet, 777.2, satisfies the impulse approximation as the timescale of the 
first stage is much shorter than the secular timescale (which is of order (777i^2/^n-+) ""^ times the 
orbital period, or ~10'^ orbital periods for jovian mass planets). Thus, the secular dynamics of the 
two surviving planets, 777-1 and 7772, would be as described in Section 2. The weak distant encoun- 
ters that eventually lead to the ejection of the third planet would produce small perturbations to 
the secular solution, as would the gravitational interactions between mi and 7773. (Neglecting the 
latter is justified, in the lowest order, as the 777.1,7773 separation is much greater than the 7772,7773 
separation.) The efi'ects of these perturbations on the secular dynamics presented above will be 
considered in a future study. 

The case of v Andromedae The analysis given in the previous section can be used to estimate the 
pre- impulse orbital eccentricity of at least one of the planets. A good candidate amongst the known 
exo-planetary systems to apply this theory is the v Andromedae system, which represents perhaps 
the most (dynamically) constrained of the known exo-planetary systems. A recent orbital solution 
for this system (as quoted in Chiang et al. 2001) gives present eccentricities of planets C and D of 
0.25 and 0.34, respectively, and yields secular dynamics for these two planets similar to that shown 
in Fig. l(e,f), with apsidal libration amplitude of 25°-35° (for relative orbital inclinations < 20°). 
An independent orbital solution obtained previously by Stepinski et al. 2000 (from an older and 
therefore smaller set of observational data) has orbital parameters within ^ 2a uncertainties of the 
more recent solution; it yields secular dynamics similar to that shown in Fig. l(a,b), with apsidal 
libration amplitude of 80°-90°. The libration period for this system is ~ 7 x 10^ yr. Now, if we 
allow that planet D suffered an eccentricity impulse, cd ~ 0.35, in its history, we can estimate 
that the initial eccentricity of planet C could have been as small as ~ 0.07 if the apsidal libration 
amplitude is 25-35 deg, but even smaller, 0-0.01, if the amplitude is 80-90°. 

After submitting this paper, we became aware of a recent preprint by Chiang &; Murray 
(2002) who propose an adiabatic eccentricity perturbation to explain the apsidal resonance in the 
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V Andromedae planetary system. They invoke torques from an exterior massive primordial disk to 

provide the adiabatic eccentricity excitation to the outermost planet D. However the adiabaticity 
of the perturbation is owed to a specific choice of disk parameters, which are not well constrained; a 
different, and perhaps equally plausible, choice of disk parameters may instead provide an impulse 
perturbation (E. Chiang 2002, personal communication); then the secular dynamics described here 
would apply. 

As mentioned previously, the apsidal libration amplitude provides an observational diagnos- 
tic: large amplitudes favor the impulse mechanism whereas small amplitudes favor the adiabatic 
mechanism. In the case of v Andromedae, there remain significant uncertainties in the orbital pa- 
rameters, and at present it appears difficult to discriminate between the adiabatic and the impulse 
mechanisms. The 25-35 deg apsidal libration amplitude found by Chiang et al. (2001) is neither 
clearly small nor clearly large; in both the adiabatic mechanism and the impulse mechanism, an 
initial eccentricity near 0.06-0.07 of planet C can account for this libration amplitude. On the other 
hand, the large libration amplitude of 80-90 deg found by Stepinski et al. (2000) would be explained 
by initial ec ^ 0.01 in the impulse limit, but cannot be explained naturally in the adiabatic limit. 
We urge further observations and analysis to improve the accuracy and fidelity of orbital solutions 
of exo-planetary systems; this would help to constrain their dynamical history. 

We acknowledge comments from E. Chiang, discussions with F. Rasio, and research support 
from NASA. 
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Fig. 1. — Secular dynamics for ji = 0.5, a = 0.33 and the following initial conditions: (a)-(c) 
eio = 0,zj7io = 0; (d)-(f) eio = 0.1e2f,ti7io = 0; (g)-(i) eio = 0.1e2f,n7io = vr. The upper two 
panels show the time variation of the eccentricities ei, 62, and the difference of apsidal longitudes, 
Atu = vji — ■uj2\ the unit of time is 27: /gi. The lower two panels plot the cumulative probability 
density function of |Ati7| and of ei. 
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(a) 




Fig. 2. — For mi = 0, the secular variation of mi's eccentricity vector is the vector sum of a 
forced component of magnitude eif along the direction of m2's apsidal longitude and a rotating 
free component of magnitude ei^free (a). For initial eio < s, apsidal resonance occurs for initial 
conditions {kio, hio) in the area of intersection of the two circles shown in (b). 
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Fig. 3. — Apsidal libration amplitude as a function of initial wio for fx = 0.5, a = 0.33, and three 
different values of eio/e2f = 0.03,0.1,0.3. 



